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ABSTRACT 

Theoretical modeling of the daytime Venus ionosphere is used to 
augment the measurements of Mariner V made during its 1967 flyby of 
Venus. The equations of heat transport for the electron, ion, and 
neutral gases are simultaneously solved along with the associated mo- 
mentum and chemical equations for the ion and neutral gas densities. 
When solutions are constrained by the observations from Mariner V, 
the resulting temperature and density profiles can be used tc estimate 
the possible magnetic field strength within the Venus ionosphere and 
concentration of neutral helium at 100 km altitude. Several alternate 
models are presented to account for the effects of eddy diffusion in 
the neutral atmosphere and solar wind heating. In all these models 
the resulting solutions lead to an approximately horizontal magnetic 
field in the topside ionosphere of strength 10 - 20y and a neutral 
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helium density of order 10 cm at 100 kin altitude. 


INTRODUCTION 


During the 1967 flyby mission, Mariner V obtained a single cross 
sectional profile of the day and nighttime Venus ionosphere (Fjeldbo and 
Eshleman, 1369; Mariner Stanford Group, 1967). This data combined with 
in situ measurements of the nightside by Venera IV can be used to con- 
struct additional features of the Venus ionosphere from theoretical 
models. The main features observed during the Mariner V encounter are 
summarized in Figure 1. At the time of observation the solar wind speed 

_3 

was about 590 km/sec with a density of 3 protons cm" and a kinetic proton 

5 

temperature of 3x10 °K. The solar wind carried a small magnetic field 

of strength about 10y oriented at 26° to the flow streamlines at the 
orbit of Venus (based on a spiral field model). Upon encountering the 
highly conducting Venus ionosphere currents are induced to flow so that 
a magnetic field appears across the front side of the planet. This mag- 
netic field deflects the inoming solar wind particles and guides them 
around to the back forming a tail -like structure extending many Venus 
radii (Dessler, 1968; Johnson, 1968). The Mariner V occultation data 
clearly shows these features in that a sharp cutoff in the electron density 
(ionopause) occurred at about 500 km on the dayside and 3500 km on the 
nightside. As shown in the figure, the location of these points was 
approximately 45° below the subsolar wind point and 135° above, respectively. 
Waves propagating upstream from the ionopause region give rise to a bow 
shock as observed in both the density and magnetic field data at about 
50,000 km from Venus along the Mariner V entry trajectory. 
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A model of the solar wind interaction with Venus and its conse- 
quences for the ionosphere has been recently discussed (Bauer, et. al . , 
1970). In this model the streaming energy of the solar wind is con- 
verted to magnetic and thermal energy in the region inside the bow shock. 
The magnetic field strength is expected to rise from about lOy in the 
solar wind to about 40y in the vicinity of the ionopause. Between the 

bow shock and the ionopause the proton temperature increases to nearly 
6 

4x10 °K. Within the ionosphere the inward solar wind pressure must be 
balanced by the sum of the ionospheric charged particle pressure and 
possible magnetic field pressure. In the ideal case of an infinitely 
conducting ionosphere, the induced currents corresponding to the solar 
wind flow lie entirely in the ionopause surface and in the region above. 
The resulting magnetic field is excluded from the ionosphere. Even with 
a realistic finite conductivity the magnetic diffusion time t^ is so much 
greater than the reversal time for the interplanetary magnetic field, 
that magnetic diffusion into a quiet ionosphere should be negligible. 
Since t^ is also much larger than a Venus day, the rotation effects and 
field reversal would have to be considered simultaneously to properly 
assess the degree of magnetic field penetration. However, if a magnetic 
field is present within the ionosphere, then it may arise from a weak 
planetary magnetic moment, from currents generated by large scale inhomo- 
genities, or from enhanced magnetic diffusion. If the effective distance 
over which the magnetic field can be considered homogeneous is reduced 
from the order of a planetary radius to a much smaller value (thereby 
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reducing t^), magnetic diffusion would become important. This may be 
possible if the likely turbulence in the topside ionosphere produces a 
wave structure with small characteristic lengths. The need for some 
source of magnetic field within the topside ionosphere is borne out by 
the model calculation presented in the following sections. 
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BASIS OF THE THEORETICAL MODELS 

From the proceeding discussion it is clear that the sun, in addition 
to the usual EUV radiation effects, strongly influences the Venus iono- 
sphere. In order to assess the effects of solar control on the Venus 
ionosphere it was necessary to solve the heat transport equations for the 
electron, ion, and neutral gases simultaneously with the momentum and 
chemical equations for their temperatures and densities. The equations 
and parameters used for this study are fully listed in the Appendix, and 
a detailed discussion of a similar set of equations as applied to the 
earth's ionosphere has been given previously (Herman and Chandra, 1969). 

In addition to the pressure balance boundary condition needed at 
the ionopause, the values of the neutral densities at the lower boundary, 

100 km, are among the most important parameters needed to characterize the 
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ionosphere. In this connection we find that a C0 2 density of 5x10 cm 
with a small amount of N 2 (<0. 1%) in photochemical equilibrium with their 
ions reproduces the lower daytime ionosphere quite well. Varying the CO,, 
density produces pronounced shifts in the altitude of the electron density 
maximum as shown in Table 1. That is, we find the calculated electron 
density maximum occurs at altitudes ranging from 120 km to 145 km as the C0 2 
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density is varied from 1x10 cm . The shift is in direct response to 
changes in altitude of the maximum solar EUV absorption. Since the top- 
side thermal structure is relatively insensitive to the C0 2 concentration 
in the range indicated in Table 1, the observed electron density peak of 

5 _ 3 

5.3x10 cm at 140 km can be used as an indirect boundary condition that 
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fixes the density of CO^ at 100 km at 5x10 cm 

The concentrations of hydrogen and deuterium at the 100 km base 
are determined by solving for all the temperature and density profiles 
simultaneously, and varying the base values of hydrogen and deuterium un- 
til the resulting densities are consistent with those inferred from 
Lyman a measurements at great distances from the planet (Barth, et. al., 

1968; Kurt, et . al , 1969). As is described elsewhere (McElroy and Hunten, 
1969), one can obtain in this manner hydrogen and deuterium densities of 
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5x10 and 2x10 cm” at 100 km, respectively. Molecular hydrogen, while 
present in the lower ionosphere, is not thought to be an important con- 
stituent for the altitude range of interest because of its rapid photodissoc- 
iation (Dalgarno and Allison, 1969). There is no experimental evidence 
concerning the quantity of helium in the atmosphere, although an estimate 
has been made based on an earthlike model of its radiogenic origin (Knudsen 
and Anderson, 1969). Accordingly, we have treated the neutral helium 
density at 100 km as a parameter of the problem, and have investigated its 
effect on the charged particle density and temperature distributions. 

The boundary conditions for the heat transport equations are speci- 
fied in terms of temperature gradients for each constituent at the upper 
boundary (500 km at the dayside Mariner V occultation point) and a common 
fixed temperature at the lower boundary (100 km). The use of a non-zero 
temperature gradient at the upper boundary implies a knowledge of the 
thermal structure above the boundary altitude. Based on a hydromagnetic 
analysis of the solar wind interaction with the Venus ionosphere 
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(Spreiter, 1970), the maximum amount of thermal energy available for 
conduction into the ionosphere is less than 10% of the energy deposited 
by energetic photo-electrons in the topside. At most this represents a 
small correction towards higher temperatures and pressures. Accordingly, 
we have used a zero temperature gradient at the upper boundary, z u , and 
have introduced a distributed source of thermal energy in the topside to 
simulate possible solar wind heating, Q CU( = (aV/A)P exp ((z - z n )/A). 
That is, a fraction a of the solar wind energy flux VP is deposited 
at each altitude z with an e - folding distance A. We have also in- 
cluded the small amount of heating from the direct coulomb interaction 
between ambient charged particles and solar wind protons that are able 
to penetrate the ionopause region (see Q g p, Q^, Q^p in the Appendix). 
Even with the maximum possible penetration of solar wind protons, the 
resulting ionospheric heating is only a 1% perturbation. 

Helium ions in the topside ionosphere are produced mainly by photo- 
ionization and lost by charge exchange with H, D, and electron re- 
combination. For rates appropriate to Venus, the daytime electron densi- 
ty, as measured by Mariner V can be supplied in about 1000 seconds, and 
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so maintained throughout the Venus day of about 10 seconds unless the 
upward velocity of He + exceeds 0.5 km/sec. A similar calculation for 
H + including both photoionization and charge exchange yields a filling 

5 

time of 10 seconds and an upper limit on the upward velocity of 0.01 
km/sec. Since the Venus day is much longer than both filling times, H + , 
D + , and He + are possible constituents of the topside ionosphere with He + , 
possibly the dominant ion. This concision depends upon the amount of 
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neutral hydrogen and deuterium present in the topside. That is, H den- 
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sities in excess of 10 cm could lead to H as the dominant ion via 
charge exchange with He + and photoionization. The only requirements 
placed on the model by the data are that the ion scale height be quite 

large (> 100 km) and that the electron density near the ionopause be 
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1.2x10 cm . This can be met by any combination of H , D , and He , 

<* _ 3 

the sum of whose densities must be 1.2x10 cm . 

An upward velocity of 0.5 km/sec supplies only 1% of the necessary 
upward pressure needed to support the topcide ionosphere against the 
solar wind. Equivalently, a supersonic velocity of about 10 km/sec is 
needed if the pressure balance is shared equally between kinetic and stream- 
ing pressure. Such supersonic velocities are not likely to be obtained 
in a system with horizontal magnetic field lines at the ionopause pre- 
venting radial outflow from the planet. When the effects of even moder- 
ate upward flow are incorporated into our models, any horizontal compon- 
ents of the magnetic field are carried with the flow and compressed into 
the topside region just below the ionopause (Banks and Axford, 1970). In 
this case we find that the electron and ion thermal conduction increases 
in the field free region thereby causing the temperatures and internal 
kinetic pressure to fall to somewhat less than half that needed for 
pressure balance with the solar wind, Figure 2. The pressure deficit can 
be made up by adjusting the parameters aV and X in Q to supply sufficient 
additional heating. These effects are discussed later. If the upward 
flow is limited by the need to maintain the observed electron density in 
the topside ionosphere, the radial density distribution of light ions is 
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approximately given by a simple force balance with V^= 0 (see 
Appendix). The solutions we obtain in this manner yield a large 
electric field in the region where the ionosphere changes from 
C0 2 + to one of the light ions, Figure 3 (E=20 m H+ e/g, where g * the 
acceleration of gravity). Newly created hydrogen ions would be 
accelerated by this large field creating an upward flux that might 
contribute significantly to the pressure balance. The result is 
that hydrogen ions would steam upwards at 6 to 8 km/sec with a flux 
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of about 10 cm sec . The flux is estimated by setting div(nV)~PL 
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where the production rate P< 1 cm sec" sec and L-100 km). The 

resulting upward pressure from such a flux f the ionopause region 

_1 1 _2 

is a maximum of lxl 0” dynes cm" as compared to the solar wind 
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pressure at the occultation point of 8x10 dynes cm . Ur ''ess there 
is a more copious source of H + in the vicinity of 200 km altitude 
than our computations indicate, it is the internal kinetic pressure 
and possible magnetic fields that supply the necessary pressure. 
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D I SCUSS ION OF MODELS: Q =0 

sw 

In support of this view, we obtained a series of temperature and 
density profiles for both the neutral and charged particle gases. The 
neutral helium density at 100 km and the effective magnetic dip angles 
were varied as parameters while all the remaining boundary conditions and 
solar EUV flux were held constant from model to model. From sets of 

4 

profiles where the electron density N g = 1.2x10 , the charged particle 
pressure, P c = N g kT g + [H + ] kT^+ + [D + ] kT D+ + [He + ] kT He+ , at the upper 
boundary was obtained and is shown in Figure 4. The solid curves repre- 
sent the kir.etic pressure, P , and the dasher’ lines the difference between 

2 

the solar wind pressure of the horizontal component, /8tt, of any mag- 
netic field that may be present in the Venus ionosphere. At low helium 
densities the pressure decreases slowly with increasing amounts of helium 

8 _ 3 

forming a pressure minimum in the vicinity of [He]^ 00 = 5x10 cm . 

Further increases cause the pressure to rise very sharply. The reasons 
for this behavior are illustrated in Figure 5 where the temperature and 
pressure curves are given simultaneously. The topside density corres- 

4 ^ 3 

ponding to each point on the pressure curve was 1.2x10 cm' . It can be 
seen that the minimum in the pressure curve arises from decreasing ion 
temperatures as the ion thermal coupling to the neutral gas increases. 
Further increases in the neutral helium density cause the solar EUV 
heating of the electrons to increase at a much faster rate than the 
corresponding loss to the cooler ions and neutral gases. Therefore, at 
low neutral helium densities the pressure follows the ion temperature 
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while at higher densities it follows the electron temperature. 

2 

The solutions appropriate to Venus are where P = P - !• /8tt. 

v# bW L 

Clearly effective magnetic dip angles less than 4° are not possible since 
the resulting electron temperature and therefore P are too high for any 

V 

concentration of neutral helium considered. Solutions for effective dip 
angles greater than 5° are not explicitly ruled out except that they re- 
quire very high neutral helium concentrations or stronger magnetic fields 
than was measured on the nightside by Venera IV (Dolginov et. al . , 1968). 

If Venus has only a weak intrinsic magnetic field (less than 10y above 
100 km altitude) that is essentially horizontal in the vicinity of the 
occuitation point, or if there exists an externally generated field of 
similar nature, then the neutral helium concentration at 100 km is likely 
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to be in the vicinity of 5x10 cm" (typical of the curve labelled dip = 4.1°). 
The total atmospheric content of He associated with this number density 

9 

can be accumulated during 4.5x10 years, as long as the total escape rate 
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of He is less than 10 cm sec. In the absence of any experimental evi- 
dence for an induced field, the Venera IV measurements of about 10y form 
an upper limit to any intrinsic field. The case B = 0 is not allowed since 
a pressure balance cannot be achieved (it is equivalent to the Dip = 90° 
case). The dashed line labelled B = 0 should be interpreted as 
0.01<B<0.05y, a field strength sufficient to reduce thermal conduction 
across field lines, yet small enough not to significantly affect the 
pressure balance (Leontovich, 1965). 

As each point on. the P c curves in Figure 4 corresponds to a separate 
set of temperature and density profiles, we show a typical example 
of one of the allowed solutions, Figure 6, for a neutral helium density 
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of 5x10 cm , P c = 8.65x10" dynes cm" , and solar EUV flux scaled from 
the values given by Hinterregger et. al . , (1963). In this model t K e 
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electron density has a peak of 5.2x10 cm at 140 km altitude in a 
region that is mainly CC^*. At 260 km the scale height rapidly increases 
as the ionic composition changes from to the light ions He + and H + . 

The profiles for H + and D + have been combined by summing. At 500 km, 
the ionopause is indicated schematically by the dashed line P = P . 

A very similar set of profiles is obtained when the quantity of neutral 
hydrogen and deuterium is increased so that He + becomes the minor ion. 

The only significant change is a further increase in scale height corres- 
ponding to the reduced average mass of the resulting ionosphere. These 
density profiles shown were computed self consistently with the temper- 
atures shown in the other half of Figure 6. Here the temperatures just 
below the ionopause are T^ = 640 °K, T g = 3670 °K T^ g + = 2160 °K, and 
T h + = 2330 °K. Because of the efficient transfer of thermal energy between 
the ions, T^+ and T He + are closely coupled in the topside. In the lower 
region where the neutral helium density is much larger, the heat loss from 
He + to He predominates over the heat input from the electrons via the large 
momentum transfer cross section corresponding to the symmetric charge 
transfer reaction. Since the coupling of H + to the neutral atmosphere 
is not as efficient as He + , T^+ tends to follow T . 

An interesting feature of the Venus thermal structure is the rela- 
tive stability of the neutral gas temperature with respect to changes in 
neutral gas composition. An example of this is shown in Table 1 where 
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the CO^ density is varied from 1x10 cm to 1x10 cm with little 
change in T . Unlike the similar case for the earth's ionosphere 
(Herman and Chandra, 1969) a single constituent, CC^, is responsible 
for both the conversion of solar radiation into thermal energy and its 
subsequent reradiation in the infrared. Increasing the concentration of 
CO2 for a given solar energy input redistributes both the energy absorp- 
tion and reradiation to somewhat higher altitudes without altering the 
exospheric thermal balance. 
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SOLAR WIND HEATING: Q ♦ 0 

The piccure of the dayside ionosphere presented so far is one where 
a relatively static pressure balance is achieved as long as the solar 
wind streaming pressure remains constant. A decrease in the streaming 
pressure would require the ionopause boundary to move out a considerable 
distance from the planet to reestablish the necessary pressure balance. 
Since the solar wind streaming pressure is known to change quite fre- 
quently relative to the flux of solar radiation, the ionopause level 
would be constantly moving. This would mean the Mariner V data was 
transient in nature, and as such is merely a snapshot of a highly dy- 
namic system. However, if the solar wind is responsible for about half 
of the charged particle heating within the ionosphere, for example by 
collisionless damping of fast mode hydromagneti c waves (Barnes, 1966) 
originating in the ionopause region, then reducing the solar wind 
streaming pressure leads to a corresponding reduction in ionospheric 
topside kinetic pressure. Using the expression Q = ( oVP /A) exp 
((z-z u )A) to simulate solar wind heating as a function of the parameters 
aV and A, the resulting electron temperatures are shown in Figure 7 and 

3 

the corresponding kinetic pressures in Table 2. When aV = 10 , the 
results are approximately the same as the case where there is no solar 
wind heating (aV=0). For example, if the heating mechanism were to 
arise from waves propagating at approximately the acoustic speed of 

5 

about 5x10 cm/sec, then the overall heating efficiency a would have to 
be 0.001 for the effect to be negligible. The overall heating efficiency 
is that fraction of the solar wind streaming energy that propagates down- 
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vertical extent of the ionosphere, A -100 km. As aV increases, the heat- 

4 

ing effect becomes very large, causing electron temperatures of 10 °K 
6 

when cA/ = 10 . From Table 2, however, it can be seen that a V cannot exceed 

4 

about 5x10 without the internal kinetic pressure exceeding the solar wind 
pressure. The values in Table 2 and the profiles in Figure 7 have been 
obtained for that density of neutral helium that produces a pressure 

4 

minimum, see Figure 4, dip = 5°. When aV = 10 , Q sw is approximately 
half the total heat input to the electrons in the topside of the ionosphere. 

7 

If V is about the ion-acoustic speed and A ~10 cm, then an efficiency of 
a = 1 to 5% would lead to a nearly stable ionopause and topside iono- 
sphere. 

In tne complete absence of a magnetic field in the topside ionosphere 
solar wind heating is insufficient to raise the electron and ion tempera- 
tures enough to maintain the ionopause. When the extreme upper limit for 

5 

solar wind heating is used, aV = 10 and A =100 km, the resulting charged 

_ 9 2 

particle pressure, P c , is only 4x10 dynes cm , a factor of two too 

4 _ 1 

small. With more likely values for Q , dl ~ 10 cm sec" , there must be 
an essentially horizontal magnetic field throughout a significant portion 
of the topside ionosphere. As mentioned earlier, moderate upward veloci- 
ties will compress any horizontal components of a magnetic field into a 
restricted region below the ionopause. If the extent of this region, d^, 
is about 100 km or larger, then sufficient heating occurs and a solution 

can be obtained. For a field of strength about 10y and dip angle of 5°, 

** _J 

oW = 10 cm sec and A= 100 km, a suitable solution is found with neutral 
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helium densities at 100 km altitude in the vicinity of 10 cm for 100 
km < d b < 400 km. 

At this time it is not possible to use the existing data to choose 
between the several models since each one leads to similar temperatures 
and densities for reasonable values of their parameters. The most likely 
choice should include a combination of all the effects we have discussed. 


t 
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THE EFFECTS OF EDDY DIFFUSION 

All of the models shown so far were computed for the case where the 

6 _2 _ 1 

eddy diffusion coefficient was taken to be 10 cm sec . Eddy diffusion 
and thermal escape was incorporated into our computations of the neutral 
atmosphere in a manner similar to that of McElroy and Hunten (1969), as 
well as the additional effect of charged particle temperature coupling. 

An example of the dependence of eddy diffusion for neutral helium is 

5 6 7 

given in Figure 8 for four values of K .. . (0, 5x10 , 5x10 , and 5x10 
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cm sec’ ). The behavior of neutral helium is as expected; in the region 
dominated by COg. the scale height for neutral helium is approximately 
that of C 02 » while at higher altitudes the more usual scaie height domi- 
nates. In the absence of eddy diffusion, = 0, the reduced scale 

height effect persists since the effects of thermal diffusion are propor- 
tional to the strong temperature gradients in this region. Hydrogen be- 
cause of its high bulk velocity does not show substantial dependence 
on either eddy or thermal diffusion. Even in the absence of eddy diffusion 
the neutral hydrogen is strongly depleted relative to a simple hydro- 
static density distribution by the upward flux escaping the planet. Thus, 
if Venus has a relatively non-., bulent atmosphere between the cloud tops 
and 150 km, helium and deuterium are likely to be enhanced relative to 
hydrogen in both the neutral and ionic forms above 100 km altitude. 

Since neutral helium and deuterium are strongly effected by eddy 
diffusion, the thermal structure of the ionosphere must also show corres- 
ponding effects. Figure 9. Above 200 km, increasing decreases the 

amount of neutral helium and deuterium present so that the ion-neutral 
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thermal coupling is reduced and the ion temperatures rise. The effect on 
the electrons is quite different. Even though the heat loss to the ions 
and neutral gases is reduced, the heat input from the solar EUV is reduced 
even more to produce a net cooling effect. Since the electrons dominate 
the charged particle pressure, the upward pressure must also decrease. 

This means that the more turbulent the atmosphere, the more neutral helium 
will be needed to maintain the pressure balance under any fixed set of 
the other parameters. 

Since the degree of turbulence in the lower ionosphere may be 
small relative to the similar height range on earth (McElroy and Hunten, 
1969), a set of pressure curves is given for no eddy diffusion, Figure 10. 

As expected, the minimum in the P c curves is shifted towards lower helium 
densities at 100 km (which become higher helium and deuterium densities 
above 250 km). In all other aspects the character of P for no eddy diffus- 

L 
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ion is very similar to the case where K . . = 10 cm sec" . For a fixed 

3 ^ 

set of parameters, B t = 10y,oA/ = 5x10 , and a sufficient ionization rate 
to maintain the topside electron density, turbulence in the bottomside 

8 

ionosphere shifts the likely neutral helium density at 100 km from 1x10 

_ 3 8 _ 3 

cm" to 5x10 cm . Corresponding to one of the curves shown in Figure 10 
(Dip = 4.1°), the topside thermal structure is shown in Table 3. Except 
for the shift in neutral helium density necessary to maintain about the 
same topside temperatures and densities, the fessential features are un- 
changed. For the same topside pressure, P c , the ion temperatures are 
slightly decreased and the electron temperature slightly increased in the 
absence of ionospheric turbulence. Therefore, while ultimately important 


t 
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for quantitative agreement with more detailed data, eddy diffusion does 
not alter the qualitative results concerning the charged particles in the 
ionosphere even though the topside ion composition and thermal structure 
depend strongly on the degree of turbulence. 
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SUMMARY 

Theoretic.^ modeling of the daytime Venus ionosphere has been used 
in conjunction with the Mariner V data to help explain some of the ob- 
served features. The observed ionopause at 500 km altitude arises from 
the deflection of the incoming solar wind by the topside ionosphere. The 
resulting pressure balance between the solar wind and ionospheric kinetic 
and magnetic pressures leads to predictions of the thermal structure within 
the ionosphere. In order to achieve the required pressure balance a 
small essentially horizontal magnetic field must be present in the topside 
to reduce the effective thermal conductivity even when additional heating 
is supplied by the solar wind. Consistent with the thermal structure and 

4 _ 3 

the observation of about 10 electrons cm in the topside ionosphere, the 

8 _ 3 

neutral helium density is likely to be in the range of "10 cm at ICO 

km. The actual value depends upon the degree of turbulence present in the 

6 2 _ 1 

lower ionosphere. With an eddy diffusion coefficient of 10 cm sec 
we find that the nearly isothermal topside electron and ion temperatures 
are approximately 4000 °K and 1300 °K respectively at the point corres- 
ponding to the Mariner V occultation data. The magnetic field for this 
case has a strength of about 10y and a dip angle of 4.2°. Such a quanti- 
tative picture of the ionosphere is only tentative since most of the 
fundamental parameters, energy inputs and boundary conditions, are 
poorly known. In particular, a good quantitative picture would require 
a knowledge of ionospheric turbulence, the degree of solar wind heating, 
the densities of the light neutral constituents, H, D, and He, and the 
effect on the dayside ionosphere of lateral transport of H + , D + , and He^ 
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Figure Captions 

1. Illustration of the solar wind encounter with Venus showing the 
bow shock and ionopause. The plane of the Venera and Mariner V 
orbits has been rotated into the plane normal to the Venus-Earth 
line. 

2. A comparison of the effects of upward ionic bulk flow (solid lines) 
on the ionospheric thermal structure with a stationary ionosphere 
(dashed lines). 

3. The electron density and ion composition when hydrogen (plus deuterium) 
is the dominant ion. 

4. The pressure of the charged particles just below the ionopause level 

(500 km) at the Mariner V occultation point vs. the neutral helium 

density at 100 km altitude. The dashed lines represent the difference 

between the solar wind pressure and the magnetic pressure of the 

? 

indicated field, P " B. /8 tt. The cases shown are for magnetic 
w t 

2 -1 

dip angles between 1 and 10°, and K e( j. = 10 cm sec . 

5. Topside electron and ion temperatures vs. the neutral helium density 
at 100 km altitude. 

6. The ion composition and thermal structure when helium is the dominant 
ion in the topside. The parameters used are a V = cm sec \ X = 100 km, 

K eddy = 1q6 cm sec ^ ’ ^ He ^100 km = 5xl ° 7 cm ’ ^ C0 2^100 km = 5x10 cm * 

t He ^T00 km = 6xl ° 6 cm 3 ’ ^100 km = 2x1 ° 6 cm ^ and ^2^100 km = 10 cm * 

-9 -2 

The resulting topside charged particle pressure is 8.65x10 dynes cm . 
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7. The effect of a nonthermal source of heating, Q. originating in 

sw 

and above the ionopause region on the electron temperature. 

Q = (oyP,./X)exp((z-z 11 )/A) . Solid lines represent A - 300 km 
and the dashed lines X = 100 km. T n is the average of the neutral 
temperatures for all six cases shown with the spread indicated by 

O C 1 

the horizontal bars. 10 <aV<10 cm sec" where a is the heating 
efficiency and V is the energy transport velocity. 

8. The neutral helium density for values of the eddy diffusion co- 

7 2 -1 

efficient ranging from 0 to 5x10 cm sec . The dashed line repre- 
sents the C0 2 density for comparison with the neutral helium scale 
height for the case where T p = 650 °K. 

9. The electron temperature for three values of the eddy diffusion co- 
efficient. P c is in units of 10" 9 dynes cm" 2 and the listed helium 
densities are calculated at 500 km altitude with a fixed value at 
the lower boundary (100 km), see Figure 8. 

10. Charged particle pressure just below the ionopause vs. the neutral 
helium density at 100 km in the absence of eddy diffusion, K eddy = 0. 
See Figure 4. The shaded region defines the range of allowed low- 
spheric solutions for the indicated magnetic find strengths. 
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APPENDIX I 


The heat transport equations for the electrons, ions, and 
neutral gases are of the form 

. ST 

v c v y- - si [w y] - « JS + - 5 v <t j ■ v 

—3 o 

N. = density (cm ), k = Boltzmann's constant (ergs/ K) , 

J 

C y = specific heat (3/2 or 5/2), Tj = temperature of species 

j(°K) , K. = thermal conductivity, Q. = solar heat production 
J J s 

_0 1 

(ergs cm sec ), Q. = energetic particle heat production 

J P 

—3 -1 

(ergs cm sec ), L .. = coefficient of the rate of energy ex- 

change between gases, z - altitude (cm), t = time (sec). 

For the electron gas: 

T 5/2 

K e ' 1 ' 8 * 10 " 5 -THT- 


Q es - S £ N j -x *x e 




■ = Y, — - — / N (z') a dz' 

\ " cos X i s v s,\ 


2 ne 

^ep ^ E 


n InA f m. 1 

f — [ir »i<V - “WJ *± (E i' 


z) 


-3 ( T) 

N. = density of jth neutral species (cm ), a. v ' = ionization 
J J > "■ 

cross section for wavelength X (Hinteregger , et. al., 1965), 

(A) 

a ' = absorption cross section for wavelength X (Hinteregger, 

J , x 

et. al., 1965), e = heating efficiency at altitude z, 

A. 

-2 -1 

0 ^ = solar EUV flux (cm sec ) (Hinteregger, et. al., 1965), 

X “ solar zenith angle, e = 4.803 x 10 -10 (e.s.u.), = energy 

of incoming energetic ion (ergs), m.^ = mass (gm) . 

-2 -1 

0 (E^,z) = flux of incoming energetic ions(cm sec ) 


2 

\ 1 . (Xj = V7 / Vt~ e-t dt = erf (X ) 


l e 


2 -X 

.4 ' (X ) = V7“ Vx„ e e 


i e 


m E . 

X = -p- — 

m i kT 


Energy loss from electrons to neutrals 


L = 1.602 x 10" 12 n (T - T )il- 77 x 10' 
en e e n \ 


N. (1-1.21 x 10~* T ) T 
1 e e 


+ 2.9 x 10 -14 N,T 1 ^ 2 + 3 x 10 -16 N„ (1-4.8 x 10 _5 T 1/2 ) 

X 6 Z 0 


+ 5.8 x 10 14 N 0 T “ 1/2 + 1.8 x 10" 15 No(l + 4.3 x 10 -4 T ) 

^ 0 O 0 
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+ 1.72 x 10 -14 NoT _1//2 + 2.46 x 10 -17 N. T 1/ ' 2 

3 e 4 e 

+ 9.63 x 10" 16 N-U-1.35 x 10 -4 T ) T 1/2 

d e e 

+ 4.82 x 10‘ 16 N g (1-1.35 x 10" 4 T 0 ) T 0 1/2 
+ 1.37 x 10" 14 N 0 exp (2.72 x 10 -6 (T - T ) 2 

£t 6 Zi 

+ 6.43 x 10~ 3 (T - T )1 

e n | 

Energy loss from electrons to ions 

L eI = 5.18 x 10' 7 n e n 1 (T 0 - Tj)/(Aj T q 3/2 ) InA 

A 1 = 28, A 2 = 44, A 3 = 2, A 4 = 4, Ag = 1 , A g = 2 

corresponding to the ionic species Ng + , C0 2 + , Hg + , He + , H + , D + , 
respectively. 

15 ^ InA ^ 20 for the cases of interest 

E L ek < T e - V - L e„ + E L el 

k I 

Ion gases: 

Only the light ions are assumed to have temperature 
significantly different from the neutral gas temperature T . 

Kj = 7.37 x 10“ 8 Tj 5/2 A z ~ 1/2 
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Energy loss from He to neutral gases 


L 4n n 4 (T 4 


- V [ : 


1.21 x 10" 26 N. (T. + T ) 1/2 
4 4 n' 


+ 5.78 x 10" 26 N x + 6.47 x 10" 25 N 2 


+ 1.28 x 10“ 25 N, + 1.09 x 10” 25 N_ + 1.17 x lO -25 

O O 


Energy loss from H to neutral gases 


L 5n ‘ n 5 (T 5 


- v[ 4 - 


078 x 10 26 N. (T. + T ) 1//2 
o o n 


+ 2.59 x 10~ 2 ° Nj + 2.04 x 10 -26 N 2 


+ 1.41 x 10 -25 N 3 + 4.69 x 10 -26 N 4 + 1.28 x 10 -25 N g J 

S 45 - 5.35 x 10- 16 — 5 ( V >4 f 1 J72 - - S 54 

= rate of heat exchange between ions 4 (He + ) and 5 (H + ) 


£ L J* (T J - V ‘ L j„ + S j,I - L ej 


Deuterium, D , is assumed to have its temperature between and 

T e . This causes little error since T. is close to T c at all 
5 4 5 

altitudes. 

Neutral Gases: 

All the neutral gases are assumed to have a common tempera- 


ture T 
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Radiative Enery Loss from Neutral Gas 
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t 

ft 


K = K . . , the coefficient of eddy diffusion 
eddy 

(McElroy and Hunten, 1969) 

„ 3 f k T n (m s + " > 1 ) 

SJ 8^1 2,1 "s m j 

a . = - 1/4 = thermal diffusion coefficient (McElroy and 


Oj = 4.0 x 10 8 a 3 = 2 x 10~ 8 

c _ = 1.06 x 10" 8 o a = 1.06 x 10 -8 

5 6 

V . (z) = N (z ) V . (z )/N (z) 

sb s esc sb esc s 

V sb ^ z esc ^ = Jeans escape velocity for species s. 

g ( 2 ) , 877 T 6-056 X10 8 I* 

[6.056 x 10° + z J 

Z L nk < T n - T k> - L nr - Z L In ' L en 

k 1 

The momentum and continuity equations have the general form 
dV i 3V i 3p i 

m i n i Tt + V i ~Tz 57 " m i n i 6 + e i n i E “2 K ij (V i " V 
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“5? + h (V i n i ) " P i ■ L pi 


Hunten , 1969) 


c ij ■ 1/2 <°i + V 


= 3.68 x 10 


-8 


ct. = 3 x 10 
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-8 
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Because of the presence of the Venus plasmapause at an 

altitude of 500 km acting as a cap on the ionosphere, the 

vertical flow velocities are likely to be small. Specifically, 
i i 4 

for { - V j | <10 cm/sec the Venns ionosphere can be well 
represented by the static equations (V., ~ 0) for the light ions, 


-5* ~ m i n iS + e i n i E " 0 


n g = n^ (charge neutrality) 


p i ■ n i kT i 


The heavier ions, C0 o and N„ are assumed to be in photo- 

2 

chemical equilibrium throughout their region of importance. 


V + C0 2 - 1 *2 * C0 2 


k 1 - 9 x 10 


C0 2 + + e - 2 CO + 0 


kg “ 3.8 x 10 


Ng + + e - 3 N + N 


kg “ 2.9 x 10 


hv + COg - C0 2 + e 


hv + Ng - Ng + e 


Hinteregger et. al (1963) 


Reactions with the light species: 

The following table summarizes the charge transfer cross 


sections and electron recombination cross sections used in this 


work. 
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Ionization 
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Energy 

Def. 
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ev 
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0 
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Banks, 1970 


co 2 

N 2 + 
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1-17 (-9) 

Calc** 
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0 
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Biondi, 1963 
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McElroy 1969 
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Biondi 1963 


* (Hochstim, 1969) 







Table 1 

PARAMETRIC VARIATION OF 
NEUTRAL CARBON DIOXIDE 


C0 2 

T e 

V 

V 

T n 

N6 P 

z p 
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1380 
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5.21x10^ 
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1550 
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5.51x10 s 
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Table 2 

The Effect of Solar Wind Heating on Charged Particle 

Pressure P„ 
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N He 

T 

e 

T He+ 

V 

T 
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5x10 
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3310 

3360 

683 

5x10 

3800 
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3310 
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5x10 6 
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2840 

2960 
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5x1 0 

3870 

1700 

1890 

680 

5x10 

5780 

856 

955 

658 

5xl0 9 

10,600 

643 

690 

625 
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Topside Temperatures vs Neutral Helium Density 
at 100 km. =0. Dip = 4.1 
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